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Abstract We present an X-ray timing analysis of the Crab pulsar, PSR B0531+21, using 
the archival RXTE data. We have investigated the stability of the Crab pulse profile, in 
soft (2-20 keV) and hard (30-100 keV) X-ray energies, over the last decade of RXTE 
operation. The analysis includes measurement of the separation between the two pulse 
peaks; and intensity and the widths of the two peaks. We did not find any significant time 
dependency in the pulse shape. The two peaks are stable in phase, intensity and widths, 
for the last ten years. The first pulse is relatively stronger at soft X-rays. The first pulse 
peak is narrower than the second peak, in both, soft- and hard X-ray energies. Both the 
peaks show a slow rise and a steeper fall. The ratio of the pulsed photons in the two peaks 
is also constant in time. 
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1 INTRODUCTION 



The Crab pulsar (PSR B0532H-21) was discovered in 1968 (IStaeUn & Reifensteinl 1 1968b . and is the 
remnant of a supernova explosion (SN 1054). The neutron star rotates about 30 times per second. Th e 
pulsed emission is visible across the entire electromagnetic sp ectrum ( radio: Hankins & Eilekl ( l2007h . 
optical: Oosterbor ek et all (|2008|) . X-rays: iMineo et all ( l2006l) . 7-rays: iKuiper et al.1 (1200 ih . and very 
high energv: iAUu et al.l ( l2008b ). 

In all the energy bands, the Crab pulse profile exhibits a dual peak structure and the pulse mor- 
phology is known to vary as a function of photon energy. The two peaks are separated by a phase of ^ 
0.4 over all the wavelengths jPravdo. Angelini & Hardingill997HKuiper et all 120031: iRots et al.Ll2004l) . 
But arrival times of the two peaks are not aligned in phase at different wavelengths . X-ray and 7-ray 
pulse leads the radio pulse (iKuiper et al.l 120031: iRots et al.l l2004t lAbdo et al.l I2OIOI) . The phase delay 
measurements between optic al and radio o bservations show a fuzzy behaviour, with reports of opti- 
cal pulse leading ( iO osterbroek et al.l l2006l) . trailing behind (iGolden et al.i , i2000i) and being coincident 
(iRomani et al.l I2OOI ). with the radio pulse. 

The variation of Crab pulse profile as a function of photon energy has been studied umpteen times 
with significant detail. The relative intensity of the two pulse peaks move in see sa w manner with in- 
creas ing energy. The first peak is more pronounced at radio and optical wavelengths (lOosterborek et al.l 
l2008b . The second peak starts domin ating as energy increases; and eventually towers above the for- 
mer at sof t 7-rays (Mineo et al.ll2b08b . With a further increase in energy, the first peak again becomes 
dominant (IKuiper et al. , 200 ll) . 

The time variability o f the Crab pulse profile is relatively less studied (IJones. Smith &Nelsonill980l; 
ICarraminana et al" l fl99l . Nevertheless, it is equally important. It can assist our understanding of the 
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emission nature of this object. Irrespective of the physical internal process in the neutron star, phenom- 
ena that can possibly change the high energy pulse profile of a rotation powered pulsar are: (i) some 
reconfiguration of t he neutron star crus t and magnetic field structure, which can be associated with star- 
quaks and glitches jpranco et al.ll2000h and (ii) ch anges in the neu tron star spin axis with respect to the 
Une of sight, which can happen due to precession ( lMelatosll2000l) . So me variations of the gamma-ray 
pulse profile of Crab was interpreted to be due to a 1 3 year precession ^Nolan et al.[ 19931 lUlmer et"ail 
[l994j) . Pulsar glitches are expected to produce thermal X-r avs tHui & Cheng. , 2004 ) and certa i n form 
of energy release can also produce pulsed X-ray emission dTang & Cheng , 200 li Icheng et al.Lll998h . 
which can change the pulse profiles for some period after the glitches. 

In this paper, we present the first ever long term pulse shape analysis of the Crab pulsar, using 
data spread over 10 years (2001-2010). Observations made with the instruments on board the RXTE 
satellite, were used to analyze the evolution of the pulse profile in — soft (2—20 keV) and hard (30—100 
keV) X-ray energy band. In section 2, we describe the data reduction from the RXTE-PCA and HEXTE 
instruments. This section also explains the model used to quantify the observations. The results obtained 
from the detailed timing analysis are discussed in Section 3. We summarize and discuss the implications 
of the observations in Section 4. 



2 OBSERVATIONS AND ANALYSIS 

Data for the present analysis were obtained from two different detectors on board the RXTE: (i) the 
Proportional Counter Array (PCA), and (ii) the High-Energy X-ray Timing Experiment (HEXTE). The 
field of view of these non-imaging detectors contained both the Crab pulsar and the Crab nebula. The 
data used in the present analysis are in the 2-20 keV band from the PCA detectors and in the 30- 
100 keV band from HEXTE instrum ent. Crab nebula suffered significant scattering during 1995-1998 
dWong et al.L boOlt iRots et al.L |2004 . Therefore, we did not include data before this period, so as to 
avoid misleading results due to broadening of the pulses. The entire analysis was done using the FTOOL 
from the astronomy software, HEASOFT-ver 5.2. The following subsections give the details of PCA and 
HEXTE data reduction. 



2.1 Reduction of RXTE-FCA data 

The RXTE-PCA consists of an array of five collimated xenon/m ethane multi anode propo rtional 
counter units (PCU) with a total photon collection area of 6500 cm^ dJahoda et al.Lll996ll2006l) . Data 
is available in different event modes and they are binned into different energy channels. For the present 
work, we have used all the publicly available data in event mode, E_250/is_128M_0_ls. This data 
mode has a time resolution of 250 yus and energy information spread over 128 channels. Selection of 
this mode, allowed us a uniformity over the energy bins over a long timebase from 2001-2010. 

From each observation, energy resolved pulse folded light curves were generated using the FTOOL- 
FASEBIN. This tool corrects the photon arrival times and converts them to the solar system barycenter 
The arrival time of each event is then converted into the phase of rotational period of the pulsar The 
output of FASEBIN is a two dimensional histogram of photon count rate against the pulse phase and 
energy. Absolute phase was dete rmined by using the time solution avail able from the Jodrell Bank 
Crab Pulsar Monthly Ephemerij^l dLvne. Pritchard & Graham-Smithlll993h and the corresponding pul- 
sar coordinates. Thereafter, another FTOOL - FBSSUM was used, which integrates the pulse profiles over 
different energy ranges. 

A typical 2—20 keV RXTE-PCA pulse profile, binned into 128 phasebins, is shown in Figure 
1 (left Panel). The figure shows the average profile of all the observations made between January to 
June, 2009. The average profile was generated using the FTOOL-FBADD. The folded light curve has an 
asymmetric, dual-peak shape. The first pulse is stronger than the second pulse. A significant inter peak 
emission (bridge emission) is also seen. Hereafter, first peak will be referred to as "PI" and the second 
peak as "P2". 



http://www.jb.man.ac.uk/ 
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Fig. 1 A sample of 2-20 keV and 30-100 keV, pulse profiles of Crab pulsar, from RXTE- 
PCA and HEXTE observations. The solid line shows the best fit model, with a reduced of 
1.3 and 0.9, for 118 d.o.f., respectively. The observation details are mentioned in each panel. 



2.2 Reduction of RXTE-HEXTE data 

The RXTE-HEXTE consists of 2 clusters each conta ining four Nal/Csl phos wich scintillation detec- 
tors, with a total photon collection area of 1600 cm dRothschild et al.Lll998h . Each of these clusters 
rock along mutually orthogonal directions to provide background measurements. When cluster-0 is on 
target, cluster- 1 is off target and vice-versa, i.e. the HEXTE rocks between on- and off- source posi- 
tions. For the present work, data in science event mode was analyzed from all the observations made 
between 2001-2010. From the raw science event files, the source file were generated, using the FTOOL- 
HXTBACK. Energy resolved pulse profiles were obtained using the FTOOLS- FASEBIN and FBSSUM. The 
FASEBIN extracts only the on-source data from the raw file. The HEXTE data was also folded with the 
same radio ephemeris, as in case of PCA data. 

Figure 1 (right panel) shows a sample of average 30—100 keV pulse profile of Crab pulsar, from 
data obtained between January-June, 2009. As opposed to the soft X-ray pulse profile, the data above 
30 keV has low statistics. We therefore, separately, generated HEXTE pulse profiles by averaging the 
data on half yearly basis. 

2.3 Model fitting 

To quantify the analysis, we fitted a 10-parameter model to the pulse profiles. The model consisted 
of position and peak intensity of each pulse peak; 2 exponentials (1 rise and 1 decay) for each pulse 
peak; and 2 constants (one for the nebular and background emission; and the other for the inter-peak 
emission). The model was fit over the phase range 0.7 -1.7. An example of the best fit model is shown 
with a solid line, in Figure 1, for both - PCA and HEXTE profiles. The reduced of the best fit ranged 
between 0.6 and 3.5, for 118 degrees of freedom. 

3 RESULTS 

All the pulse profiles of RXTE-PCA and HEXTE, were fitted with the model as described above. The 
variability/stabiUty of the important parameters of the fit, are discussed below. Several glitches have 
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Fig. 2 The separation between the two pulse peaks of the Crab pulse profile, from RXTE- 
PCA and HEXTE observations. The solid line in the bottom panel shows the half yearly trend 
in HEXTE. Arrows in the top panel indicate the time of occurence of glitches. 



been detected from the Crab pulsar in the period of these X-ray data. The times at which the glitches 
occured (taken from lEspinoza et all (1201 Ih -) are marked in the top panel of each subsequent figure. All 
except one glitch had similar magnitudes {jS.v/v of ^10 x 10~^), while one of them (which occured 
around MJD 53067) had a magnitude of ~200 x 10"^ 

3.1 Separation of the two pulse peaks 

As shown in the sample pulse profiles in Figure 1, two clear peaks are seen at phases ^^0.99 and ~1.39. 
Figure 2 shows the separation between the two peaks, from the RXTE-PCA and HEXTE observations. 
The X-axis in this figure, gives the MJD of respective observation, and the Y axis gives the phase 
separation between the two peaks. 

A constant fit to the RXTE-PCA curve gave a phase separation of 0.4058 (5), with a reduced 
of 0.94 for 189 d.o.f. A constant fit to the HEXTE results, gave a phase separation of 0.4079 (6), 
between the two pulse peaks. The reduced of fit was 0.36 for 188 d.o.f. The solid line in the results 
from HEXTE detector, show the half yearly trend. This indicate that the pulse peak maximas are stable 
in time. The separation of the maxima of the two peaks is constant over the last ten years. Further, 
the phase separation betw een the two peaks do not show a significant variation w ith energy, consistent 
with earlier observations (IWills et al.lll982l:IClear et al.lll98'il:lNolan et anil993h . We did not find any 
considerable effect of occurence of glitches on the separation between the pulse peak maximas. 

3.2 Pulse widths 

The rise times and the decay times of PI and P2 (in phase units), for the PCA and HEXTE pulse profiles 
are shown in Figure 3, along with occurence of glitches in the top panel. The X-axis gives the MJD of 
each observation. In case of HEXTE, the half yearly average trend is shown with a solid line. A constant 
was fit and results are given in Table 1, with 90% confidence limit, and of fitting. In both the energy 
bands, the two pulses are asymmetric. Both the peaks, PI and P2, show a slower rise and a steep fall. 
The first pulse peak is narrower than the second peak. It rises in almost half the time taken by the second 
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Fig. 3 The rise and decay times of the two peaks of the Crab pulse profile. The left panels 
show the results from RXTE-PCA observations, while the HEXTE results are shown in the 
right panels. The half yearly trend of the HEXTE observations is shown with a solid line in 
the figure. The arrows indicate the time of occurence of glitches. 



peak. The first peak also decays faster than the second peak. We did not find any detectable effect of 
occurence of glitches on the pulse widths. 



3.3 Pulsed intensity 



Figure 4 shows the ratio of the integrated pulsed photons in the two peaks, for RXTE-PCA and HEXTE 
observations. A constant fit to the two ratios, gave a value of 0.633 (1) for PCA and 0.604 (13) for 
HEXTE energy bands. This result is co nsistent with results reported using COMPTEL observations 



at me dium 7-ray energies (1-10 MeV) (iBuccheri et al.l 1993t Carraminana et al. . 19941: Kuiper et all 

3, 



I2OOII) : and from EGRET observations above 30 MeV dFierrol ll 995c .TompkinTet al.L 1 1 997h . However" 
more variation is seen in case of HEXTE data, as compared to PCA results, which might be due to 
poorer statistics in case of HEXTE data. Further, although the second peak is smaller, it contains large 
number of pulsed photons. 



3.4 Pulse Intensity Ratio 

A comparison of variation in intensity ratio of the amplitude of the two pulse peaks, in PCA and HEXTE 
profiles, is shown in Figure 5. Here, the intensity ratio of the two peaks (P1/P2) is plotted as a function 
of time. A ratio of 1 or above, indicates that the first peak is dominant, whereas a ratio of less than 
1 imply a stronger second peak. Clearly, PI is significantly stronger than P2 in soft X-ray band (PCA 
data). A constant fit gave an intensity ratio of 1.63, with a reduced of 0.8 for 189 d.o.f. The intensity 
ratio of the HEXTE pulses shows a variation, which could be due to poorer statistics. The solid line in 
the figure, shows the half yearly average trend in the intensity ratio of HEXTE profiles. A constant fit in 
case of HEXTE results gave an intensity ratio of 1.18, with a reduced of 0.33 for 188 d.o.f. Clearly, 
the first peak dominates at higher energy also, albeit the ratio is smaller 
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Fig. 4 The ratio of integrated counts in PI to P2, for RXTE-PCA and HEXTE observations, 
as a function of time. The arrows indicate the time of occurence of glitches. 
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Fig. 5 The variation of intensity ratios of the two pulse peaks (P1/P2) of Crab pulsar, are 
shown as a function of time. The curve at the top show the variation for PCA observations; 
while the curve below corresponds to HEXTE observations. The solid line in the lower curve, 
shows the half yearly average trend in the HEXTE observations. The arrows indicate the time 
of occurence of glitches. 



4 DISCUSSION 

A coherent picture of the pulsed emission from Crab pulsar has been an enigma inspite of numerous 
multi wavelength studies. The radiation in radio, optical. X-rays and high energy 7-rays is known to 
differ considerably. In this paper, we have reported that the X-ray pulse profile of the Crab pulsar is 
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Table 1 Results for model fitting of the Crab pulse profile. 



Parameter 


RXTE-PCA 


RXTE-HEXTE 




Value 


Reduced 


Value 


Reduced 






(for 189 d.o.f.) 




(for 188 d.of.) 


Phase separation between two pulse peaks 


0.4058 (5) 


0.94 


0.4079 (6) 


0.36 


Rise time of PI (in phase unit) 


0.03077 (5) 


1.3 


0.0335 (6) 


0.43 


Decay time of PI (in phase unit) 


0.02180 (4) 


3.5 


0.0236 (12) 


2.5 


Rise time of P2 (in phase unit) 


0.0642 (1) 


0.75 


0.0707 (15) 


0.44 


Decay time of P2 (in phase unit) 


0.03051 (8) 


0.52 


0.0241 (5) 


0.49 


Ratio of pulsed photons in two peaks (P1/P2) 


0.633 (1) 


0.63 


0.604 (13) 


0.39 


Pulse intensity ratio (P1/P2) 


1.6344(35) 


0.8 


1.18(2) 


0.33 



*The numbers in bracket, give the error for 90% confidence range in the respective values. 



stable in time over the last ten years of RXTE performance. There is no indication of time variability 
of the pulse shape at soft and hard X-ray energies, during the last decade of observations. No evidence 
of time variation in the phase separation between the two pulse peaks, the widths of the two peaks and 
their relative intensities, was noticed. 

The light curves of X-ray pulsars reflect the geometry of the magnetic field and the radiation emis- 
sion pattern of the rotating object. If the pulses are produced near the neutron star surface, the com- 
pactness (mass to radius ratio) of the neutron star will also cause gravitational bending of light and thus 
effect the pulse profile. There are several reports on the long-term flux variability of the Crab pulsar 
This is conjectured to be due to slow precession of the spin axis of the neutron s t ar, which results in 
precession of the pulsar magnetosphere at similar periods dOzell Il989l : iKanbachi 1 19901 : iNolan et"an , 
119931: lUlmer etalill994 . Since the X-ray and 7-ray emission is modulated by the magnetosphere ge- 
ometry, therefore, changes in the respective pulse profiles should be in sync. Flux variability can also 
be at tributed to sudden glitch-induced changes between the pulsar spin axis and the magnetic moment 
axis jLinket al.L 1992, 1997: Franco et al .. 2000). A change in the alignment angle during a glitch, can 
modulate the duration of the Une of sight's traverse through the pulse emissio n cone. This in turn can 
produce considerable differe nces in the pulsed emission from the neutron star dJones. Smith & Nelsonl 
\mO- Sekimo to etami995[) . 

However, if these theories are true, then the pulsed X-ray and the 7-ray emission are supposed to 
show similar variations. The 7-ray pulse profile of the Crab pulsar shows considerable variations over 
time. The i ntensity ratio of the two peaks in the pulsar's light curve, is known to undergo a periodic 
behaviour dWills et all 119821: lOzelL ll99U iNolan et al.L 1 19931 lUlmer et al. I [19941) at high 7-ray energies. 
However, we have found that the intensity ratio of the two peaks is stable in time, both in soft and hard 
X-rays. This is consistent with earlier works in X-ray band dWills et an,ll982l) . It has been proposed that 
the periodic variability in the intensity ratio of the two peaks is masked at lower phot on energies, due 
to onset of plasma-do minated radiation, which mainly affects the second pulse peak dHasinger et all 
11984 lHasing"e3.ll985h . But neither the magnitude of the radiation component is certain, nor are there 
supportive conclusions from the optical and radio observations. 

The emission of thermal radiation from the neutron star is also affected by the turmoil in its interior. 
It could be due to the rotational dynamics of the superfluid in the neutron star's inner crust. The X-ray 
light curve can a lso get modulated du e to release of large amount of X-ray energy during the sudden 
onset of glitches dHui & Chengll2004l) . However, we have seen that the pulse profile is constant in time. 
There are no significant changes in the pulse shape parameters during glitches. Therefore, it is difficult 
to comment on the extent of turbulence caused by the glitches. 

The underlying physics involved in the emission mechanism of Crab pulsar is not completely un- 
derstood and coordinated multi-waveband study of the pulsar is crucial. Our results fill in the gaps 
between the trend seen at other wavelengths and we hope that observations by more sensitive instru- 
ments on future missions, like ASTROS AT ( (PauL .2009.) . will provide important clues regarding the 
pulsar emission mechanism. 
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